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Phosphate minerals such as those in the apatite group tend to be the dominant forms of 
phosphorus in minerals on the Earth's surface. Phosphate can be reduced to phosphides 
during high-energy events, such as lightning and impacts. Here we show that, in addition to 
formation of metal phosphides, a new compound was formed by lightning in a fulgurite from 
New Port Richey, Florida, USA. A calcium phosphite material, ideally CaHPO3, was found in 
spherules mainly consisting of iron silicides that formed by lightning-induced fusion of sand 
around a tree root. This phosphite material bears a phosphorus oxidation state intermediate 
of that of phosphides and phosphates in a geologic sample and implicates phosphites as 
being potentially relevant to other high-energy events where phosphorus may partially 
change its redox state, and material similar to this phosphite may also be the source of 
phosphite that makes up part of the phosphorus biogeochemical cycle. 
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perceived as being limited to phosphates under terrestrial 

conditions. Cosmically, P-mineralogy includes phosphide 
minerals in which phosphorus mainly forms binary compounds 
with a metal such as Fe or Ni. The mineral schreibersite, 
(Fe,Ni)3P, is an example of such a binary solid and is cosmically 
the most abundant phosphide. Terrestrial occurrences of schrei- 
bersite and other phosphides have since been discovered within 
pyrometamorphic rocks!?, and especially within fulgurites*+, 
which are glasses formed by cloud-to-ground lightning. The 
presence of phosphides within fulgurites parallels the existence of 
phosphides within impact rocks such as tektites°, with impact 
cratering also being the likely route to forming phosphides in 
lunar rocks®, and in impact melts on other meteorites”®. 

Phosphorus has thus far been found in mineral form in only 
two oxidation states: as +5 in phosphates, and as ~—1 in phos- 
phides. This latter oxidation state is based on the binding energy 
of inner electrons of the P atom in phosphides?, which is con- 
firmed by ab initio calculations!°. Notably absent from known 
P-mineralogy are any other oxidation states between these, 
wherein the P anion is a reduced phosphorus oxyacid. Salts with 
P in a + 1 oxidation state (hypophosphite), +3 oxidation state 
(phosphite), and +4 oxidation state (hypophosphate) have been 
synthesized but are not known from naturally-occurring geologic 
samples. Despite the apparent lack of reduced phosphorus oxy- 
acids in mineral form, reduced phosphorus compounds are not 
rare in the geo-environment. For example, the gas phosphine 
(PH3) has also been recognized to be a trace gas in many 
environments!!-!3, and has even been implicated as a bio- 
signature gas on Venus!4. More specifically to the case of 
potentially mineral-forming P-compounds, between 1 and 10% of 
all microbially lineages can use phosphite (and often hypopho- 
sphite) as sole P sources!>, demonstrating these ions have sources 
in the environment!®. In addition, reduced phosphorus oxyacids 
have been found in geothermal spring water and anaerobic water 
in Florida where their origin is likely microbial!”!*, Within rock 
samples phosphite has been found in extracts of fulgurites!?, in 
serpentinites”°, and in extracts of Archean rock samples”!?, 
though the mineral source of such a phase is unclear. Therefore, 
reduced P oxyacids must exist in solid form of some sort in the 
environment, either as a pure compound such as CaHPOs, or 
with the phosphite exchanging for other ions such as phosphate 
or sulfate. The failure to identify a specific mineral form of 
reduced phosphorus oxyacid in the environment is problematic 
given the widespread evidence for phosphorus redox in 
biogeochemistry!®, but may be due instead to an incomplete 
search for such material. 

Due to their rapid formation which results in lithologies that 
may vary significantly with respect to temperature, pressure, and 
Joz over a few millimeters of distance, fulgurites are an ideal 
material in which to search for reduced P-compounds. Phos- 
phides such as schreibersite have been identified in multiple 
fulgurites: from Michigan*, Pennsylvania!?, and Illinois*. In 
addition, even more reduced material such as iron silicides are 
recognized in multiple fulgurites*?-7°. Such material is pre- 
sumed to form from the intensely reducing environment possibly 
generated by a combustion reaction of in situ organics*?°. Pasek 
and Block!” also reported phosphite in extracts of fulgurites, 
which were attributed in one fulgurite to phosphides, which react 
with water to release phosphite as a metastable solute*!, but in 
others the source of phosphite could not be wholly identified. 

Here we identify (1) a novel form of crystalline P within a 
fulgurite from New Port Richey, Florida, USA; and (2) propose 
that the reduction of phosphate to more reduced forms takes 
place through at least two separate pathways and that both may 
be relevant to the planetary distribution of phosphorus***. The 


T= mineralogy of the element phosphorus (P) is generally 


New Port Richey (hereafter NPR) fulgurite was obtained from 
private sellers who found the fulgurite on their property in a 
residential neighborhood in New Port Richey, Florida (28.248 N, 
82.718 W) the summer of 2012 after a thunderstorm. We then 
purchased the fulgurite in late 2012. The fulgurite formed in a 
partially drained sandy soil, composed primarily of quartz sand 
and clay*4, and was about 500 g in total mass. In addition, as the 
area receives heavy rain, roots are often coated with an iron oxide 
plaque?5-37, which is composed of a 1-2-mm-thick rind of iron 
oxide that cements quartz sand particles (Supplementary Fig. 1). 

The fulgurite is a type II fulgurite*> with a massive, glassy wall 
that surrounds the central void (Fig. la). The fulgurite’s color 
ranges from white to dark gray. The most unusual feature of the 
fulgurite is the presence of large (~1cm diameter) spherules 
(Fig. 1b) of a material with a dull gray metallic luster that are 
embedded in the glassy matrix of the fulgurite. These spherules 
have a density of 4.9-5 g/cm? and are not attracted to a magnet. 


Results 

The New Port Richey fulgurite bears both silicate and reduced 
(silicide) compositional regions. The silicate fraction of the New 
Port Richey fulgurite is volumetrically more common and is 
composed of glass with two major compositions (Fig. 1c): one 
exclusively SiOz, likely relict quartz sand grains, and the other 
composed mostly of Ca, Al, and Si oxides (Supplementary 
Table 1). Within the silicate lithologies are small droplets 
(<15 um) composed of P-rich iron metal (Fig. 1d). Such droplets 
have been reported previously in other fulgurites**1%. 

The large metal spherules within the NPR fulgurite consist of 
silicides, specifically with the composition of FeSi and FeSi, in 
addition to a material that is darker in contrast in back-scattered 
electron (BSE) imaging that bears Ca and P and is amorphous to 
X-ray diffraction (Fig. le, f). The silicides also bear 0.45-1.25 
wt.% P. The groundmass of the spherules is composed of FeSi, 
(84% by volume), in which euhedral crystals of FeSi (11%) grow, 
often together with the CaP material (5%). The nominal cationic 
composition of the spherule is thus Sig3Fe33Ca2:P2. The CaP 
material bears vermiform FeSi with a width of ~1 um. The CaP 
material is most closely associated with the FeSi material, which 
grows as euhedral crystals (30-70 um) that are larger than the 
average FeSi crystals in the FeSi, groundmass, or more compact 
(with less extreme maximum compared to minimum dimen- 
sions). Though the CaP material is in contact with the FeSi, 
groundmass, crystal growth of FeSi, within the CaP material is 
less common. Inside the CaP material there are very rare (less 
than 20 um sized) crystal fragments that turned out to be a new 
crystalline phosphite (see the dashed circle in Supplementary 
Fig. 2). 

Given the similarity in chemical composition and unit-cell 
values, the structure of the new phosphite was refined starting 
from the atom coordinates given by Phillips and Harrison’? for 
the P4;2,2 structural model of synthetic CaHPO3. Phillips and 
Harrison*® describe the three-dimensional structure of synthetic 
CaHPO; as having CaO, capped distorted trigonal prisms, with 
pseudopyramids of HPO; sharing corners and edges. The for- 
mula of the phosphite phase in the NPR fulgurite is (Cag.9;Feo 99) 
HPO3. Although there is a thin boundary between anthropogenic 
and natural fulgurites (see below), here we first report a phosphite 
compound which represents the reduced version of monetite 
(CaPO3OH), as the phosphite (HPO;2~) anion bears a +3 oxi- 
dation state for P. To the best of our knowledge, CaHPO; is also 
the first fulgurite-occurring crystalline material wherein P is not a 
phosphate or a phosphide. 

We attempted to replicate the CaHPO; mineral formation by 
adding CaHPO, to a mixture of Fe and Si (0.1:1:2 molar ratio) 
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Fig. 1 The New Port Richey fulgurite images and microscopy. a Glassy tubes that consist of a glassy melt surrounding an internal void, in turn surrounded 
by a crust consisting of cemented sand grains. The diameter of the fulgurite is 2 cm, and length is 7 cm. b Spherules of gray, metallic material within the 
fulgurite with diameters of 1.1cm (left) and 0.5 cm (right). c BSE image of the glass of the NPR fulgurite. Varied lithologies of the glass include a darker 
material (left) composed of SiOz, and a lighter material (Ca, Al-rich silicate). Within this glass (red rectangle is the region expanded) are d spherules of iron 
metal enriched in phosphorus. e BSE image of the large metallic spherules of the NPR fulgurite. These consist of FeSi (medium gray), FeSi (light gray), and 
a Ca-P-O material that includes CaHPO3 (dark gray). f The Ca-P-O material is mostly in contact with FeSi. 


and heating to 1000 °C for 8h. The phosphate was demonstrated 
to form reduced P-compounds upon extraction with NasEDTA 
solution by NMR (nuclear magnetic resonance) spectroscopy 
(Fig. 2) that demonstrates the formation of P?+ and P4* species. 
However, such chemistry naturally arises from the extraction of 
P-compounds from Fe3P (or schreibersite), and the ratio of these 
species (1:5:5:1 of hypophosphate: phosphate: phosphite: pyr- 
ophosphate, see ref. 3°) implicates that—rather than reducing 
CaHPO, to CaHPO3—reduction of CaHPO, to Fe3P occurred 
instead. If we assume that the reaction velocity doubles with every 
30°C increase in temperature (corresponding to an activation 
energy of 20 kJ/mol), then at 1220 °C—where FeSi, would have 
melted—we would expect to see the reduction of CaHPO, to Fe3P 
over the course of 5 min. Given that this did not occur in the NPR 
fulgurite, that implies the reduction timescale was less than 5 min. 

The NPR fulgurite bears a crystalline CaHPO; material that 
may be the first mineral bearing a phosphite (HPO, P?+) 
group to the best of our knowledge. Notably, this fulgurite is the 
second to have formed from New Port Richey, FL, USA (con- 
trasted with that reported in ref. 4°). However, fulgurites 
can occur in material that is also influenced by either (1) artificial 
target material such as pavement or conductor metal, or 


(2) formed through downed powerlines that introduce energy to 
the environment slowly and over time*?. Notably, the NPR ful- 
gurite shows no indication of the former as there is no aluminum 
or copper in a metallic material, and the groundmass of the 
fulgurite is consistent with local soil compositions*+. The source 
of the discharge to be an artificial source may be unlikely due to 
the incomplete reduction of phosphate to phosphide within the 
FeSi/FeSi, droplets, as reduction should be complete if the system 
attained equilibrium. The reduction experiments described above 
show CaHPO, to Fe3P would likely have occurred as tempera- 
tures exceeded 1220°C (the melting point of FeSi,) under rela- 
tively short timescales (~5 min) for this fulgurite, presuming a 
reaction rate doubling every increase of temperature of 30°C. 
That the reduction of CaHPO, to Fe3P did not occur implicates a 
shorter (<100s) heating timescale, consistent with a lightning 
strike. 

A metallic silicide spherule belonging to the NPR fulgurite was 
studied by high-resolution multiscale nanotomography. Fig. 3 
reports the color-coded 3D thickness distribution of the pore 
network inside the cement sample (see Supplementary Video 1) 
together with the thickness distribution plot of the analysis (inset 
in Fig. 3). The porosity also a directionality to it, which may be 
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the same as the orientation of the flow of gas out of the fulgurite. 
Figure 4 shows a virtual cut of the spherule evidencing melting 
and recrystallization structures. 

The reported presence of much larger iron silicides#!4? in 
geologic specimens suggests that large spherules are formed by 
lightning and may be the source of these silicides that have been 
reported in various locations. For these reasons, we argue the 
CaHPO, could be a natural material and therefore represent a 
member of a new mineral group, one predicted to occur in the 
environment from biology and biochemistry!>. 
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Fig. 2 31P NMR of Fe-Si-CaHPO, (1:2:0.1 by mole) mixture heated to 
1000 °C for 8h. The NMR spectrum, in which individual anions are 
referenced to a frequency of 161.907 MHz for H3PO, (O ppm), identifies P 
species based on parts per million deviations (8) from this frequency. The 
peaks correspond to (left to right) hypophosphate (P4+), phosphite (P3+), 
phosphate (P°+), and pyrophosphate (P°+). Note that phosphite is a broad 
doublet with a J-coupling constant of 570 Hz when coupled to 'H. About 
7% of the phosphorus in this NMR spectrum has been changed from 
HPO,2-, which is consistent with reduction of phosphate to phosphide. 
Phosphide, specifically Fe3P, reacts with water to form a 1:5:5:1 mixture of 
hypophosphate: phosphate: phosphite: pyrophosphate. In contrast, a 
change to CaHPO3 would yield only phosphite. This suggests that reduction 
of phosphate to phosphite occurs much faster than the timescale of 8h 
reaction, proceeding to Fe3P within 8h. 
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Discussion 

The NPR fulgurite presents an unusual mineralogical paradox: 
the more reduced phosphorus mineral (P!~ as a metal phos- 
phide) is in contact with a silicate glass, but the more oxidized 
phosphorus material—the calcium phosphite (P?+)—is in direct 
contact with extremely reduced iron silicides. Such a difference is 
especially exemplified by the fo2 mineral buffer diagram (Fig. 5) 
that demonstrates that the formation of CaHPO; should not 
occur in contact with FeSi or FeSiz. We examine the cause of this 
disparity in P chemical oxidation state and the implications of its 
formation. 

The presence of silicides within fulgurites has been reported 
previously?®. The size of the crystals of the silicides is comparable 
to the Winans Lake fulgurite’, wherein crystals of silicides 
~10-100 um wide were found within metal spherules. Composi- 
tionally, these silicides also match those reported in the Houghton 
Lake fulgurite”®, in which FeSi and FeSi, formed from a eutectic 
melt within spherules. The size of the silicide-bearing spherules is 
much larger than those found previously? by about a factor of 
10x, which were previously the largest spherules reported in the 
literature from fulgurites. The size of these spherules (1 cm), 
however, is comparable to those reported as purported distal 
impact ejecta*!. 

The formation of the iron silicide spherules in the NPR ful- 
gurite likely occurred during the combustion of roots coated with 
iron oxide plaque (see Supplementary Fig. 3 for details on the 
reduction of Fe,O3 and SiO, by organics). Iron oxide plaques 
form around roots, and cement sand as lithotrophic bacteria 
oxidize Fe(II) in wetlands environments*°, and such plaques are 
known to concentrate phosphate’, including as calcium phos- 
phates such as brushite*4. During the lightning strike, the com- 
bustion of organic root matter likely drove extreme reduction of 
the iron oxide plaque, melting it as FeSi, formed (T > 1220 °C), 
still well below the melting point of calcium phosphate 
(>1600 °C*5). Calcium phosphate, entrained within the iron sili- 
cide, was subjected to the extreme reduction, and locally formed 
the CaHPO; material. Given that the CaHPO; is in contact with 
FeSi, and FeSi is formed under somewhat more oxidizing con- 
ditions than FeSi, (Fig. 5), perhaps the reduction of CaHPO, to 


Pore Thickness Distribution 


Fig. 3 Nanotomography of the iron silicide spherule. The scan in one field of view allows to visualize the color-coded 3D thickness distribution of the pore 
network inside the sample. The thickness distribution plot of the analysis is shown in the inset. 


4 COMMUNICATIONS EARTH & ENVIRONMENT | (2023)4:70 | https://doi.org/10.1038/s43247-023-00736-2 | www.nature.com/commsenv 


COMMUNICATIONS EARTH & ENVIRONMENT | https://doi.org/10.1038/s43247-023-00736-2 


ARTICLE 


Fig. 4 Virtual cut of the iron silicide spherule. The spherule shown in Fig. 3 has evidence of melting and recrystallization structures. 
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Fig. 5 Oxygen fugacity diagram. Oxygen fugacity is buffered by minerals 
(MH magnetite/hematite, WM _ wustite/magnetite, |W iron/wustite) 
that demonstrates the reduction of calcium phosphate to calcium phosphite 
(CaHPO3 - CaHPO,) should occur at much higher foz than other reduction 
reactions, including the transformation of schreibersite (Sch, nominally 
Fe3P) to whitlockite (wht, nominally Ca3(PO,)2). The transformation of 
whitlockite to schreibersite should also occur readily if silicates are reduced 
to silicides. 


CaHPO, was accompanied by the oxidation of FeSi, to FeSi and 
SiO,. In addition, as Fe2+ is ~9% of the divalent cation of the 
CaHPO; material (Cao9;Feo99HPO3) iron as FeSi may have 
served as a reductant as well. Elemental mapping of the CaHPO3 
material (Fig. 6) demonstrates that both Si and O are present 
within this substance. Experimentally heating CaHPO, with Fe 
and Si at 1000 °C did not produce CaHPO; but instead produced 
Fe3P after 8h, suggesting that rapid heating and cooling of this 
material—such as what occurs in a lightning strike on the time- 
scale of 60 seconds or less—is necessary to make CaHPO3. 

The NPR fulgurite bears two distinct, reduced P phases- the 
CaHPO; in the FeSi/FeSi, spherule, and P dissolved within Fe as 
much smaller droplets within the silicate glass. These likely 
represent two separate reduction pathways for phosphorus. The 
P-rich iron spherules in silicatic glass (Fig. 1c, d) are smaller 
(15 um) and their spherical shape and higher P concentration 
along the edges of the spheres*® indicate that molten iron reacted 
with P at high temperature (T > 1538 °C), probably as volatilized 
P, through phosphidation of the metal. The spherical shape 
(Fig. 1d) also indicates that iron melted. Metal phosphidation has 
been proposed to be a kinetically favored route to phosphate 


reduction that may be the source of metal phosphides on the 
Moon’, This route involves volatilization of P from phosphates 
to form P and P, gases under reducing (but not extremely 
reducing) conditions, followed by reaction of P gas with metal to 
form metal phosphides, a process that occurs rapidly*®. 

In contrast, the CaHPO,; mineral likely formed from the 
reduction of calcium phosphate, possibly as either amorphous 
calcium phosphate, or brushite (CaHPO,:2H,O). Brushite can 
dehydrate to monetite (CaHPO,) with heating and does so 
rapidly (timescale of <1 min for similar Mg-phosphates*?) at 
temperatures >200 °C4>50, Brushite reduction and dehydration, 
as: 


CaHPO, : 2H,O = CaHPO, + 2H,O + 1/20, 


was approximated using limited thermodynamic data!? but 
should occur under more oxidizing conditions than at the 
wustite-magnetite fo2 buffer. In contrast, the FeSi or FeSi, oxi- 
dation buffer? is well below the iron—wustite buffer, just slightly 
above the Si-SiO, buffer. This suggests that the reduction of 
phosphate to phosphide proceeds through a phosphite inter- 
mediate, which did not complete due to the relatively short 
duration of heating (<1 min) during the lightning strike fusion of 
the soil. The CaHPO; compound is thus out of equilibrium with 
its surrounding assemblage and should have been reduced to 
phosphide. 

To this end, the reduction of phosphate proceeds via two 
routes (Fig. 7): a high temperature (>1700 °C) reduction wherein 
Fe metal forms and is rapidly phosphidized by P gas at moder- 
ately low fo2, and a lower temperature (>1220°C) reduction 
wherein P as calcium phosphates is entrained within extremely 
reduced material and is rapidly reduced to phosphites but kine- 
tically inhibited from forming the ultimate reduced materials, the 
phosphides. 

Although the production of CaHPO; from CaHPO; in contact 
with FeSi likely has few direct counterparts in natural geologic 
samples (as silicides are rare), the identification of a likely natural, 
crystalline form of a P+ material is important to understanding the 
role of redox in the phosphorus cycle. Hess et al.* specifically call 
out the lightning reduction of phosphate to have been a widespread 
phenomenon on the early Earth potentially transforming thousands 
of kg of phosphate to phosphides and phosphites annually, which 
would have rivaled meteoritic sources. These materials may have 
been prebiotically important*}>!. Pasek and Block!® argue lightning 
would have been a less important (~3000 kg/yr) source of reduced P 
on the modern Earth, but the relatively facile reduction of phos- 
phate to phosphite (Fig. 5) indicates that phosphites may be present 
in other geologic or cosmic samples subjected to reduction but that 


COMMUNICATIONS EARTH & ENVIRONMENT | (2023)4:70 | https://doi.org/10.1038/s43247-023-00736-2 | www.nature.com/commsenv 5 


ARTICLE 


COMMUNICATIONS EARTH & ENVIRONMENT | https://doi.org/10.1038/s43247-023-00736-2 


Fig. 6 Elemental maps of the Ca-P-O material. The Ca-P-O material (BSE image in a) is shown in three-element RGB maps, with (b) Fe as red, O as 
green, and Si as blue, and (c) Fe as red, O as green, and P as blue. The Ca-P-O material is associated with some Si and O, suggesting its formation may be 


coupled to FeSi or FeSi oxidation. 
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Fig. 7 Schematic of the two routes to phosphate reduction within the NPR 
fulgurite. The high-temperature reduction is rapid due to gas-solid 
reactions, and occurs closer to the plasma discharge, and the low- 
temperature reduction occurs in the solid phase and is likely kinetically 
inhibited from proceeding to form phosphides. Known temperature points 
are the melting point of the Ca-Al-Si matrix (1195 °C), the melting point of 
FeSiz (1220 °C), the melting point of iron metal (1538 °C), the melting point 
calcium phosphate (>1600 °C), and the melting point of SiO2 (1710 °C). 
SiO, is glass near the interior of the fulgurite, whereas it is a crystalline 
material on the exterior, showing a radial temperature profile that has been 
reported previously2°. 


do not reach very high temperatures. For instance, phosphates 
within pallasites or iron meteorites may equilibrate with iron metal 
to transform into phosphite minerals. Phosphite is also found 
within serpentinites*°, and phosphites such as CaHPO; (specifically 
the Mg/Fe?+ equivalents) may be the mineral carrier of phosphite 
that yield phosphite during the extraction of these rocks. The 
chemistry of P within reduced igneous rocks (at or below the WM 
buffer) may include reduced P material such as CaHPO3;, but few 
speciation-specific studies have been done of P chemistry in these 
types of rocks. 

Modern metabolic pathways implicate Pt as being a non- 
negligible portion of the total phosphorus inventory, but the 
absence of solid P3+ materials as sources of this ion for life has 


drawn questions as to the environmental phosphite reservoir!>. 
Thus a careful investigation of the anionic character of P within 
reduced material—though not necessarily as reduced as investi- 
gated here—might reveal other forms of crystalline reduced P 
minerals, many of which may have been important in the 
development of life on the Earth°?, or in the supply of nutrients to 
modern environments!®. 


Methods 

SEM. The instrument used was a Zeiss-EVO MA15 Scanning Electron Microscope 
coupled with an Oxford INCA250 energy-dispersive spectrometer, operating at 
25 kV accelerating potential, 500 pA probe current, 2500 cps as average count rate 
on the whole spectrum, and a counting time of 500 s. Samples were sputter-coated 
with 30-nm-thick carbon film. X-ray maps were collected with acquisition times of 
10 ms per pixel. 


X-ray diffraction. Studies were performed at the CRIST, Centro di Studi per la 
Cristallografia Strutturale, Department of Chemistry, Universita di Firenze, Italy. A 
small fragment (size about 14 x 13 x 11 um) was extracted from the polished sec- 
tion (Supplementary Fig. 2) under a reflected light microscope and mounted on a 
5- um-diameter carbon fiber, which was, in turn, attached to a glass rod. Single- 
crystal X-ray diffraction data were collected by means of a Bruker D8 Venture 
equipped with a Photon II CCD detector, using graphite-monochromatized MoKa 
radiation (A = 0.71073 A), exposure time of 30s per frame and a detector-to- 
sample distance of 7 cm. Intensity data were integrated and corrected for Lorentz- 
polarization and absorption with the Bruker software packages. A full Ewald sphere 
was collected up to 20 = 70° at room temperature. 


LA-ICPMS. Laser Ablation Inductively Coupled Plasma Mass Spectrometry was 
performed using a New Wave Research (esi) laser (NWR-213) attached to a Perkin 
Elmer Nexion 2000P Quadrupole ICPMS. Two samples of groundmass glass were 
analyzed by LA-ICPMS. Glass samples were chosen as fragments of the NPR 
fulgurite, then sequentially ground and polished to a 0.05-um alumina grit and 
were mounted using putty. Four points on each glass were ablated and measured 
for the abundance of 2Si, 48Ti, 27Al, 23Na, 57Fe, 55Mn, #3Ca, 39K, 31P, and 1B 
(each chosen to minimize the effect of other isobars), which were referenced to 
USGS glasses BCR-2G, BHVO-2G, and BIR-1G*? mounted in epoxy to constrain 
major element abundances that were normalized to 100% by weight. 


Thermodynamic calculations. Data for the Gibbs free energy of various species that 
determine oxygen fugacity as buffers (Fe.O3, Fe;O4, FeO, Fe, FeSi, FeSi,, CaSiO3, 
Ca3(PO4)2, Fe3P, SiO2, CsH}20¢, CO2(g), and O2(g)) were retrieved from HSC 
Chemistry (Outokompu Research Oy, v. 7.1, see details in refs. °4°°), Data for the 
CaHPO;/CaHPO, are from Pasek and Block!? and were calculated by determining 
the difference in solubility between CaHPO; and CaHPO, (monetite), then assuming 
the Cps of both were similar, allowing extension of the buffer to high temperature. 
While such a method is likely to be less accurate at high temperature, it provides an 
approximation of the foz buffer that may be accurate enough on the logarithmic scale 
to make conclusions about the relative importance of various reactions. 
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CaHPO, reduction experiments. CaHPO, (99% from Fisher Scientific, 0.031 g) 
was mixed with Fe and Si powder (0.1 g each) to give a 0.1:1:2 mixture of CaH- 
PO,:Fe:Si by mole. This mixture was heated in an alumina crucible under N; gas at 
1000 °C for 8 h. The resulting powder was extracted with a solution of Na,EDTA in 
water (10 mL) for 1 week, and analyzed by 3!P NMR (details in4°) for 10,000 scans 
in 'H-coupled mode. 


Nanotomography. High-resolution multiscale nanotomography was carried out by 
means of a Bruker Skyscan 2214. The sample was fixed on the sample holder with 
dental wax, and the entire volume was scanned in one field of view to visualize 
pores and different density phases inside the sample. The Paganin Phase Retrieval 
algorithm was applied to the projection images to enhance the different density 
feature inside the sample. 


Data availability 

Data generated or analyzed during this study are included in this published article (and 
its supplementary information files; https://www.researchgate.net/project/New-Port- 
Richey-Fulgurite). Raw data for charts (NMR fid, Laser Ablation ICPMS data, oxygen 
fugacity vs. temperature data for buffers, and the thermal maximum temperature profile) 
may be found publically at https://www.researchgate.net/project/New-Port-Richey- 
Fulgurite. 
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